Impact of turbulence on the stratified flow around small particles 



o 

(N 



Itzhak Fouxon 1 and Alexander Leshansky 2 ^ 

1 Department of complex systems, Weizmann Institute of Science, Rehovot 76100, Israel and 
2 Department of Chemical Engineering, Technion, Haifa 32000, Israel 

We study the turbulent flow of the density-stratified fluid around a small translating particle. It 
was found recently [A. M. Ardekani and R. Stocker, Phys. Rev. Lett. 105, 084502 (2010)] that 
without turbulence, the familiar Stokes flow is dramatically altered by the stratification. The latter 
"turns on" the buoyancy introducing a new "cutoff" length scale for the flow around the particle, 
yielding closed streamlines and a faster decay of velocity. This result, however, didn't account for 
the potential role of background turbulence, intrinsically present in many aquatic environments, 
mixing the density and thus opposing the above effect. Here we derive and solve the advection- 
diffusion equation that describes the interplay of turbulent mixing, diffusion of the stratifying agent 
and buoyancy. We derive an exact expression for fluctuations due to weak background turbulence 
and show that stronger turbulence can completely change the flow around the particle. Therefore, 
the account of background turbulence is obligatory in a typical marine environment. 



(N 

i 

q 

O 

^ ■ 

Or 



> 

00 

in 

m 

o 



PACS numbers: 47.27.-i, 47.55.Hd, 47.63.mf, 83.10.-y 

Recently it was found that density stratification can 
change the incompressible flow around the small parti- 
cles strongly In spite of wide separation of scales of 
stratification (kilometers) and swimmers (0.1 — 1 mm) 
in typical marine environment, one cannot neglect the 
spatial inhomogeneity due to stratification. The reason 
is that inhomogeneity opens a new pathway of interac- 
tion - the buoyancy. The latter has no influence on the 
incompressible flow when the distribution of the strati- 
fying agent is considered uniform, but it does when in- 
homogeneity is taken into account. A combination of 
buoyancy, diffusion (conduction) and viscosity creates a 
new length-scale L of the order of 1 mm. Beyond that 
scale the perturbation flow around the particle decays 
faster than the Stokes flow of the unstratified flow that 
holds at scales much smaller than L. In sharp contrast 
with the unstratified case, the flow possesses toroidal ed- 
dies and closed streamlines at scales of order L. It was 
suggested that this may affect propulsion of small organ- 
isms and sinking of marine snow particles, diminish the 
effectiveness of mechanosensing in the ocean |l| , stifle nu- 
trient uptake of small motile organisms [2|_or potentially 
hinder the drift-induced biogenic mixing [3J. 

Clearly, to apply this result one has to determine its 
domain of validity. However, this effect of stratification 
should be very sensitive to the presence of a background 
flow that would redistribute the stratifying agent, alter- 
ing the buoyancy forces and the resulting flow. In ap- 
plications this demands, firstly, to include into consid- 
eration the background turbulent flow that is present in 
natural environments invariably. Turbulence mixes the 
fluid opposing the effect of stratification. Like in the 
case of stratification, it is incorrect to use the separation 
of scales of turbulence (the Kolmogorov's scale £ v , see 
4],) and of L as an argument to neglect flow inhomogene- 
ity. Turbulence at scales smaller than £ v is a large-scale 
chaotic flow. In contrast to mixing in the inertial range, 



turbulence mixes at small scales at a rate independent of 
the scale 0-0) ■ That rate is given by the characteristic 
value A of the gradient of turbulent velocity field. Note 
that A 2 is the energy dissipation per unit mass e divided 
by the kinematic viscosity v. The mixing occurs both 
due to turbulence and diffusion/conduction. Turbulence 
dominates the mixing at scales larger than the Batchelor 
scale £d (where the diffusive time-scale is of order A -1 ) 
and at smaller scales the diffusion/conduction prevails. 
Thus if turbulence is strong, £d < L, then the flow de- 
viates strongly from the one in the quiescent fluid. In 
contrast, for weak turbulence £d ^> L the stratified flow 
[H holds at the relevant scales smaller or of order L. 

This Letter describes consistent, quantitative analysis 
of the flow around small particles in the presence of both 
the density stratification and the turbulence. The trans- 
lating particle distorts the surrounding turbulent flow. 
This distortion redistributes the stratifying agent, such 
as temperature, which feeds back the flow through buoy- 
ancy. We derive the advection-diffusion equation that 
describes this interaction by a consistent reduction of the 
full system of hydrodynamic equations. The main modifi- 
cation relative to the usual equation on the passive scalar 
field mixed by turbulence Q is the non-trivial wave- 
number dependence of the diffusion coefficient. Solving 
the equation we describe the flow and temperature fields 
around the moving particle. The solution depends on 
one dimensionless parameter j3 = L 2 /£ 2 d . Turbulence is 
negligible at f3 <C 1, but at (3 1 the streamlines cor- 
responding to the Stokes' flow without stratification are 
recovered. It should be emphasized that the results are 
obtained without modeling the statistics of turbulence 
and they can be applied to natural environments directly. 
In the latter, /3 > 1, where the effect of turbulence is of 
order one, is typical, so turbulence is important. Thus 
velocity inhomogeneities at large scales influence the flow 
around small-scale objects. The relevant parameter is not 
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the spatial scale £ v , but the temporal scale A -1 . 

We use the Boussinesq approximation [9( to describe 
the interaction of the flow v with the stratified agent 9 , 

d t v + v ■ Wv = -V (p/p) + 0g + vV 2 v, V • v = 0, 
d t 9 + v ■ V0 = kV 2 (9; v(\x-Y\ = a,t) = V (I) 

where pj p is the pressure divided by the density, g = —gz 
is the gravitational acceleration, Y, V and a are the 
particle's coordinate, velocity and radius respectively 
[8J). We decompose the flow into the background tur- 
bulent flow u, P, 9q and the perturbation induced by 
the boundary condition describing the particle, v(x) = 
u(x) + w' [x - Y(t)\, p/p(x) = P (x) + P [x - Y(t)\ and 
9(x) = 9 Q (x) + [x - Y(t)}, where w', P, decay at 
large r = x — Y(t). The perturbations obey, cf. Q, 

d t w' + w' ■ Vw' + [u(r + Y [<]) - V(t)] • Viu' + aw' 

= -vp + ®g + isv 2 w', d t ® + w' ■ ve + w' ■ w 

+ [u (r + Y [t]) - V(t)] ■ V6 = kV 2 0, V • w' = 0, (2) 

where we introduced ov, (i) = VjUi[Y(t), t). The bound- 
ary conditions are decay far from the particle and 
■u/(|r| = a) = — (u— V). We make several assump- 
tions. The experimentally relevant situation corresponds 
to L <C i v , so we assume the latter inequality and study 
the scales r <C £ v (this implies a « I,). Since turbulent 
velocity is smooth at scales smaller than £ v then we 
have u (r + Y [t]) mu(Y [t]) + ar. Thus u (r + Y [t]) - 
V(t) pa ar - U, where U = V(t) - u(Y(t),t) is the 
particle's velocity relative to the local flow. We assume 
Ut v /u < 1 and Ui^jn < 1 so the terms including U can 
be dropped from the equations at the scales of interest 
r <C £ n - As Xi^/v ~ 1 holds for turbulence, the first 
assumption means U is not much larger than the char- 
acteristic velocity u v ~ \l n of the viscous scale eddies of 
turbulence Q. The second assumption gives U < u^/Pr, 
where Pr = v/n is the Prandtl number. As it will be- 
come clear below at Pr < 1 turbulence is irrelevant at 
L <C Irf- Below we assume Pr 1 / 2 ^> 1, as it also fits many 
situations in aquatic environments. Then U < u v /Pr 
implies U <C u^. The latter inequality means that dur- 
ing the correlation time A -1 of a the particle's deviation 
from the trajectory of the fluid particle is much smaller 
than £ n . Therefore the statistics of a(t) is the same as 
the one of the velocity gradient of turbulence in the fluid 
particle's (Lagrangian) frame. 

The above assumptions imply that i?e rc i ~ Ua/u and 
Ua/n are small, which allows to drop all terms quadratic 
in the perturbation in the equations. Further, we as- 
sume that in the domain of interest the correction to the 
gradients of the stratified agent due to the turbulence is 
negligible and one can approximate V#o by a constant 
—jz like in the fluid at rest (note that we use the units 
where 9 is dimensionless so 7 = (l/po)(— dp/dz) has the 
dimensions of inverse length). Finally we notice that at 



r <C £ v all terms in the LHS of the first of Eqs. ([2]) can 
be dropped in comparison with the viscous term. In con- 
trast, the terms dt® and ar ■ V8, that are both of the 
order A8, are comparable with kV 2 at the diffusive 
scale id = (ft/A) 1 / 2 = i^Pr^ 1 ^ 2 -C i v and are not neces- 
sarily negligible. Thus <9 t 9 + (ar ■ V)0 should be kept 
at r <C £ v . Summarizing the above, we obtain 

VP = ®g + vV 2 w', d t O + (ar ■ V)8 - ~fw' z = kV 2 9. 

Finally, to study the flow at r > L where L ^> a, one can 
model the effect of the translating particle on the flow by 
a point-force term in the momentum equation, 



VP = <dg + vV 2 w' + fS(r)z, V ■ w' = 0, 
d t Q + (ar ■ V)6 - jw' z = kV 2 6, 



(3) 



where I is the unit vector in the upward direction [15| . 

Turbulence is described by dt® + (ar ■ V)0 term which 
is the usual term describing the advection of the passive 
scalar fields by turbulence at large Prandtl numbers, see 
7] and references therein. This term is comparable with 
the diffusive term at a characteristic scale i& — 
and it is dominating at larger scales. This is the term 
which effect we study in present work. 

We briefly consider the case L <C id where to leading 
order the advection term can be dropped at r < L. This 
is the case without turbulence considered in Taking 
the Laplacian of the first of Eqs. © we find, 
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-(gj/K)w' z + vV 4 w' + V 2 fS(r). (4) 



The squared Laplacian term and the buoyancy term 
g~fw' z /K are of the same order at the scale L = (isn/jg) 1 / 4 
introduced in At smaller scales viscosity dominates 
and the Stokes flow holds. At r ~ L a non-trivial change 
in the flow pattern around the particles occurs [H. 

We return to the full system (0}. Taking the Fourier 
transform and using the incompressibility tr a = yields 

ikP = &g- vk 2 w' + fz, k ■ w' = 0, (5) 
d t Q- (a l k- V)0- 7 u/ = -nk 2 ®. (6) 

We multiply Eq. ([5]) with k and use the incompressibility 
condition ik ■ w' — to eliminate the pressure, 



P = igk z ®'/k 2 , ®' = ®-f/g. 



(7) 



Introducing k = k/k and the projection Iljj(fe) leads to 

vk 2 w' = ®'U(k)g, = (8) 

We note that obtaining Eq. ([7|) involves division by k 2 
making it necessary to consider the point k = sepa- 
rately, see below. Substituting w' z in Eq. ([6]), we obtain 
the following closed advection-diffusion equation for 

d t ® - (a l k ■ V)0 = -a(k)Q + (j>(k), a(k) = nk 2 d(k) 



(9) 
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where k\ — k 2 — k 2 . The stratification produces wave- 
number dependent diffusion coefficient d(k) and the 
source of the fluctuations (f>. Although both quantities di- 
verge at k = 0, the time evolution of O(fc) at k ^ decou- 
ples from k = 0, which is clear from relations below, so 
one can consider the solution at k and then continue 
it. Note that Eq. (|9|) has the solution 0(fc) = constxS (k) 
at zero stratification 7 = that describes constant dis- 
tribution of in the real space. One finds that 0' obeys 

d t e' - (a*k ■ V)0' = -a{k)e' - Kk 2 f/g. (10) 

To find the solution we pass to the moving frame 
Q(k,t) = Q'(k(t),t) where k(t) = W-^\t)k with 



W = aW, W 



-l.t 



W ij (t = 0) = S ij . (11) 



Since a is statistically the same as the velocity gradient 
of u in the fluid particle's frame, then W is statistically 
the same as the Jacobi matrix of the turbulent flow back- 
ward in time That is, if we consider the Lagrangian 
trajectories q(t, r) defined by dtq(t, r) = u[t, q(t, r)} and 
q(t = 0,r) = r, then Wij(t, r) = djqi(t,r) at t < de- 
scribes the evolution of small volumes in the turbulent 
flow backward in time and obeys Eq. ([lip. In particular, 
since the Lyapunov exponents of the backward in time 
flow are (—A3, — A 2 , — Ai) where (Ai, A 2 , A 3 ) are the Lya- 
punov exponents of the forward in time flow, then k(t), 
which is governed by W~ l,t rather than W(t), obeys 



lim (l/\t\)ln[k(t)/k(0)} = Ai, 



(12) 



see details in [7]. Thus the growth of k(t) with \t\ is 
similar to the exponential growth of the separation be- 
tween two infinitesimally close fluid particles in turbu- 
lence (governed by the principal Lyapunov exponent Ai ) . 

The limit in Eq. (fT2|) holds for almost every realization 
of a(t) and does not involve the randomness of turbulence 
that disappears after taking the infinite time limit. To 
describe the fluctuations of k{t) when t is finite, one in- 
troduces the polar representation k(t) = k exp[p(t)]n(t) , 
where |n| = 1. Using k = —o t k one finds Q 



-cr*n + p 



c, c 



(13) 



It follows that \n[k(t)/k] = /" ((t')dt' where C is a finite- 
correlated noise which correlation time r c is of order of 
the correlation time of a, so that r c ~ A" 1 . Thus Eq. (fT2|) 
resembles the law of large numbers. To find the moments 
of k(t) one introduces 



lim (l/|i|)ln<fc'(t)> 

t— y— 00 



4,(1). 



(14) 



The function (f)(1) is convex. It obeys (f)(0) = 4>(— 3) = 0, 
so it is negative at —3 < n < and positive otherwise. 
This holds independently of the statistics of turbulence 



(see Q for details). In the moving frame Eq. (|T0| be- 
comes 



d t e = -a[k(t)]Q - Kk 2 (t)f/g. 



(15) 



We consider at t = 0, taking the initial condition at 
t = —T and studying the limit T — > 00, i. e. we focus on 
the steady state solution. Using Q(t = 0) = Q'(t = 0), 



e'=- 



dt cxp 



a[k(t')]d£ 



k 2 (t). (16) 



The above together with Eqs. ©-© give implicit solu- 
tion to the system (j3|) in the Fourier space. Since for 
turbulence k(t) is a random vector, then and w are 
random too and should be studied statistically. The com- 
putation of the statistics however cannot be done due to 
the complex dependence on a. Thus we consider the lim- 
iting cases, which will allow us to understand the behav- 
ior of the solution in detail. We introduce new integration 
variable Sk(t) = f. a [k(t')] dt' that depends on k via the 
final condition fe(0) = k. We have 



Q'(k)=-Lr 

gJo 1 



exp [— s] ds 



■k 2 ± (s)/L 4 k e (s)' 



(17) 



where k(s) = k[tk(s)]. We introduce the scaling form 
tk/h(s) = T k (ps)/X, where d s r k = -k 4 (t j '\)[k\(T / 'A) + 
^(t/X)]- 1 and t(s = 0) = 0, so that 

,0 

s = / k 2 (t'/\)[l + kl(t'/X)/k 6 (t , /X)]dt'. (18) 
Introducing q(s) = fc[rfc(s)/A] one finds 



f f°° 

- / F(Ps)e- s ds, F = 
9 Jo 



ql+q & 



This form of the solution is particularly well-suited for 
the study of the impact of turbulence because F(s) is 
determined by turbulence only. If intermittency (depen- 
dence of the statistics of a/X on the Reynolds number 
Re) can be neglected, then F(s) varies at scales of order 
one. The solution depends on one dimensionless parame- 
ter, j3. Studying how it changes when (3 is increased from 
(the stratified flow without turbulence) to 00 one finds 
the impact of turbulence. Intermittency is negligible up 
to rather high Reynolds numbers of 10 5 — 10 6 due to the 
smallness of the corresponding anomalous exponents [4] , 
so we consider this case first. Thus F((3s) in Eq. (|T9|) 
varies at scale 1//3. The limit of small f3 is described by 

00 

0' (k/L) = -(f/g) P n F^(0), = 0' + fig. (20) 

n=0 

To find the solution to order /? one can use —F'(Q) = 
F 2 (0)(1/F)'(0) = F 2 (0)[q 2 ± /q 6 }'(0). Writing the result 
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in terms of the solution without turbulence Oo (k/L) 
(f z /g)k\ [k\ + fc 6 ] _1 fl, one finds 



We can use 



-2p(t)\ 



Q(k/L) = Q (k/L) 



1 + SQ (k/L) 



(21) 



SO (k/L) = 2p\- 1 k & kJ 2 [fc 6 +kj] 2 f ■ a f k, 

where / = [(fc 2 - 3k 2 ) k x , (k 2 - 3k\) k y , -3k 2 ± k z ] . 
Thus in this order the relative correction 5Q to the solu- 
tion without turbulence is a linear function of the current 
value of the gradient of the turbulent velocity V jUi at the 
location of the particle. Since the statistics of a is close 
to the Lagrangian statistics of V jUi (the statistics in the 
frame of fluid particle), then (<Jij) = (Vjiti) = and 
(cij<7 mn ) = (VjUiV n u m ), where the last average can be 
taken in Eulcrian frame due to incompressibility, giving 



30v(aija mn ) 



SijSmn] • (22) 



where isotropy (typically valid for small-scale turbulent 
fluctuations that determine Vit Q) and spatial unifor- 
mity are assumed. The form of (VjWiV„u m ) is fixed 
uniquely by the demands of isotropy, incompressibil- 
ity and spatial uniformity that implies (V m UiV„u m ) = 
(V m [uiS/nUm)) = 0. The relation is exact due to sta- 
tionarity condition v '(V ju/SI jUi) — e. Thus we obtain 
the exact result for fluctuations of 9 (clearly (SO) = 0) 
around 0o that are caused by the turbulent fluctuations 
of the background velocity field 



( ( 58 2 ) 1 /2 =/3fc 8 fc -2 [ fc 6 +fc 2J V[4fc2/2 - 2(k • /) 2 ]/15. 



The result indicates that the impact of turbulence is of 
order one at j3 ~ 1. When intermittency is taken into 
account, the higher order terms in the series (|20|) . that 
holds in the limit (3 — > irrespective of the neglect of 
intermittency, have weight increased by a power of Re. 
Thus the impact of turbulence is of order one at /3 ~ 
Re" 5 , where the phenomenological exponent 5 is small. 

We now consider the limit of strong turbulence (3 = 
\(v / n^g) 1 / 2 3> 1 when intermittency is negligible or Re 
is fixed. One observes that this limit is equivalent to 
the one of small stratification 7 —> in agreement with 
consideration that strong mixing opposes stratification 
canceling its effects completely at (3 — > 00. We consider 



gQ(k/L) 
f 



qj_(/3s)e- s ds 



qi(s)( 



qi(f3 S )+qS(f3s) J Q (3[ql(s) + q^s)} 

The RHS is equal to one at fc = recovering 8(0) = / / g 
following directly from Eq. ([5]) . If fc ^ then the leading 
order behavior is obtained by setting exp[— s/ f3] = 1, 

gO(k/L) _ 1 f°° ql(s)ds 
f PJa qi(s) + q%s) 

Returning to the original integration variable, t — 

Tk(s)/\, 



o(l/ /3). (23) 



q 2 ± (s)ds 



ql(s)+q 6 (s) 



A f° e - 2 ^[n 2 x (t) + n 2 y (t)]dt. 



exp[0(— 2)t] following from 
Eq. (fT4l) at Xt >• 1. Since n(t) and p(t) can be consid- 
ered as independent, while n is distributed isotropically 
0, we find (e~ 2p ^n 2 (t)) = exp[0(-2)t]/3. Thus 

(e(fc/L))~2Aj7[3/3 3 fc 2 <X-2)]. (24) 

The integral at fc = is convergent due to exp[— s//3] 
while the integral at fc ~ 1 is obtained by setting 
exp[— s//3] 1. It can be seen that the transition be- 
tween the two asymptotic regions is where the two an- 
swers are of the same order, that is (3k 2 ~ 1. The result 
in terms of O' is particularly simple, 

&'(k)^-L, fc 2 »-^; l©'(fc)l<-> fc 2 «i- 



g 



L 2 p' 



g 



L 2 (3 



For finding the inverse Fourier transform, one can put 
0' « — f/g at large (3 uniformly. Substituting 0' = — f/g 
in Eq. (|5]) recovers the Stokes flow vk 2 u)' = fU(k)z. 

Thus the solution depends on one dimensionless pa- 
rameter f3 so that the Stokes flow with streamlines open 
everywhere holds at /3 3> 1 and the stratified flow with 
closed streamlines holds at j3 <C 1. It follows that at 
(3 > 1 turbulence cannot be disregarded. In particular, 
there is a critical value f3 c ~ 1 such that the streamlines 
open at (3 = f3 c (there are closed streamlines aty (3 < f3 c 
but not at f3 > j3 c . The account of intermittency (the 
dependence on the Reynolds number) is expected to pro- 
duce stronger fluctuations of a and 0. It should lower 
the value of (3 where the account of turbulence is neces- 
sary so that in the limit of very high Re turbulence can 
be important already at /3 <C 1. 

Thus turbulence is important at j3 > 1. It can also be 
important at /3 <C 1- In this case turbulence is negligible 
in w' , but the total flow v = u + w' can differ from w' 
significantly due to u. The flow u produces the charac- 
teristic difference XL of velocities of particles separated 
by L. The corresponding relative velocity induced by w' 
is estimated by Ua/L. The ratio of the two differences 
XL 2 /Ua = (v/Ua)(L 2 /t 2 ) is the product of the large pa- 
rameter v/Ua = l/Re re \ and the small parameter L 2 /£ 2 . 
If it is large turbulence is important. Thus turbulence is 
important at (v /Ua)(L 2 / 'f 2 ) > 1 or |/3| > 1. 

Let us now estimate the typical values of (3 in various 
aquatic environments. Using the extreme value of the 
density gradient 7po = 1 kg m -4 [l[ that may occur lo- 
cally in fjords [T3|, lakes and reservoirs [l3| with \i = 1CP 3 



0.6 mm for salt-stratified water 
2 mm for temperature- 
1.4 x 10~ 7 m 2 s _1 ). Further con- 



9 m 2 s : ) and L 



kg m s yields L 
(k « 1.3 x 10 
stratified water (k 
sidering weakly turbulent conditions with the dissipa- 
tion rate per unit mass e « 10~ 10 m 2 s -3 (e.g. Kunzc 
et al.[ll| measured e < 10~ 9 in a coastal inlet) gives 
A = \Je/v ss 0.01 s _1 . Thus, the corresponding values 
of f3 ~ XL 2 1 k are « 0.3 and ss 2.8 for temperature- and 
salt-stratified water, respectively. Furthermore, in the 
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marine environment the buoyancy frequency N = ^/gj 
corresponding to the marginal oscillations which the sta- 
ble stratification supports [3] is typically in the range 
between 10~ 4 and 10 -2 s _1 , yielding density gradients 
7Po that ranges between 1CT 6 and 10~ 2 kg m -4 , several 
orders of magnitude lower than that considered in [l| . In 
some extreme cases, however (e.g. during seasonal thcr- 
mocline [lH) N may exceed 0.05 s^ 1 so 7^0 may reach 
~ 0.3 kg m~ 4 . Using this extreme value of density strat- 
ification and e sa 10 -10 m 2 s~ 3 we arrive at j3 f=s 0.5 
and P ~ 5.5 for temperature- and salt-stratified water, 
respectively. However for the less extreme conditions of 
marine turbulence and/or stratification typically j3 > 1. 
For example, for e 10~ 9 m 2 s~ 3 and jpo ~ 0.01 kg 
m~ 4 we find j3 ss 8 and f3 w 90 for temperature and salt 
stratification, respectively. 

We derived and solved the advection-diffusion equation 
that describes the turbulent flow around small translat- 
ing particles in the stratified fluid in the limit of large 
Prandtl numbers. We showed that when intermittency 
is negligible (which allows very high Re) the solution is 
determined by one dimensionless parameter /3, where tur- 
bulence is important when ft > 1. Intermittency, impor- 
tant at higher Re, only strengthens the impact of tur- 
bulence. We conclude that the account of turbulence is 
necessary in natural environments. 
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